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Inflationary models can correlate small-scale density perturbations with the long-wavelength 
gravitational waves (GW) in the form of the Tensor-Scalar-Scalar (TSS) bispectrum. This correlation 
affects the mass-distribution in the Universe and leads to the off-diagonal correlations of the density 
held modes in the form of the quadrupole anisotropy. Interestingly, this effect survives even after 
the tensor mode decays when it re-enters the horizon, known as the fossil effect. As a result, the 
off-diagonal correlation function between different Fourier modes of the density fluctuations can be 
thought as a way to probe the large-scale GW and the mechanism of inflation behind the fossil effect. 

Models of single held slow roll inflation generically predict a very small quadrupole anisotropy in TSS 
while in models of multiple fields inflation this effect can be observable. Therefore this large scale 
quadrupole anisotropy can be thought as a spectroscopy for different inflationary models. In addition, 
in models of anisotropic inflation there exists quadrupole anisotropy in curvature perturbation 
power spectrum. Here we consider TSS in models of anisotropic inflation and show that the shape 
of quadrupole anisotropy is different than in single held models. In addition in these models the 
quadrupole anisotropy is projected into the preferred direction and its amplitude is proportional 
to QtNe where Ne is the number of e-folds and g, is the amplitude of quadrupole anisotropy in 
curvature perturbation power spectrum. We use this correlation function to estimate the large scale 
GW as well as the preferred direction and discuss the detectability of the signal in the galaxy surveys 
like Euclid and 21 cm surveys. 


I. INTRODUCTION 

The recent cosmological observations such as Planck 2015 are in very good agreement with the inflationary cosmology 
as the source of the density perturbations in the CMB and the origins of structures in the Universe [T]. Although CMB 
observations are consistent with the single field slow-roll (SFSR) inflation models but there can be other types of the 
fields including additional scalar, vector and tensor modes which could play a role during inflation in a distinguishable 
way [5]. Therefore, it is worth to search for their observational imprints. Among the above additional fields, tensor 
modes are very interesting because they are one of the inevitable predictions of general relativity. There have been some 
experiments searching for the imprints of the gravitational waves in the CMB [Si! however there was no detection of 
primordial tensor mode. From the Planck results an upper bound r < 0.1 for the ratio of the amplitude of tensor 
perturbation to scalar perturbation in CMB power spectrum is obtained. While the CMB experiments are near to 
their completion, it is worth to search for primordial GW in large-scale structure surveys like Euclid and 21-cm surveys. 
There were also some theoretical works looking for the imprints of the GW in the GMB and large-scale surveys such as 
lensing by the tensor modes of the galaxy distribution, 21-cm fluctuations etc. ISHIH. 

The key point in the galaxy searches is the two-point correlation function or power-spectrum. Under the null 
hypothesis, this two-point function is statistically homogeneous and isotropic. In addition, there is no correlation 
between different Fourier modes. Coupling the inflaton field with the tensor mode leads to off-diagonal correlation 
between different Fourier modes and thus can be thought as a way to probe the primordial GW. However, since the 
final goal of the galaxy surveys is to look at the power-spectrum of the density perturbation, it is convenient to assume 
that there is a hierarchy between the scale of the density perturbations, supposed to have small wavelength, and that 
of the tensor mode, assuming to have a long wavelength. Therefore one should consider the three-point function of the 
Tensor-Scalar-Scalar (TSS) at the squeezed limit. It turns out that this three-point function looks like a quadrupole 
anisotropy in power spectrum. 

The next step is considering the non-linear evolution of this bispectrum during the radiation and matter dominated 
eras and finally, in order to read the observed mass distribution, it is also required to consider the projection effects 
P51 - E5] . Intuitively, this effect comes from the fact that in a redshift survey one can find the position of the galaxies 
by using their apparent position and redshift in the sky. However, in the presence of the tensor mode, the apparent 
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position will not be the same as the original position of the object. As a result there are three different potentially 
important contributions that must be taken into account in order to evaluate the observed power spectrum of the mass 
distribution in the Universe. However, among the above effects, only the first one depends on the specific inflationary 
model. So one can neglect the non-linear and projection effects and consider only the primordial bispectrum. 

It has been shown that for SFSR inflation models, all of the above effects are at the same order and there would 
be a cancellation between them. This seems to be the case for all models which satisfy the Maldacena’s consistency 
condition [55]. As a result, for these kinds of models, the signal would be small and beyond the accessibility of the 
current galaxy surveys. So it is worth to go beyond the SFSR inflation and see how the violation of the consistency 
relation shows up in the observed power spectrum. This has been done very recently in [27H2Hj . As they have shown 
there are no cancellation between different components in these models. The authors have claimed that the quadrupole 
amplitude is observable by the current galaxy surveys. However, its shape is exactly the same as in SFSR inflation 
models. 

It is well-known that the GW has a nontrivial evolution, an oscillatory and decaying profile, when it re-enters the 
horizon. However, it can be shown that the above phenomenon is a cumulative effect and does survive even after the 
long wavelength tensor mode decays inside the horizon [^ . As a result this effect is called the fossil effect. 

In this paper, we generalize these studies to models of anisotropic inflation in which there is a preferred direction 
during inflation yielding to statistical anisotropies in the form of quadrupole anisotropy in curvature perturbation 
power spectrum. In addition, the bispectrum has a different shape and the TSS would be aligned along that direction 
in contrast to SFSR models. In models of anisotropic inflation we assume that at the background level there is a vector 
field which sources the background anisotropy. As we shall see in the following this anisotropy affects the TSS and we 
will obtain 


(^k'l^CkaCka) = (x • e(^^(ki) • X - kg • XX • (ki) • kg^ (ka)(ki) + (kg ka) (1) 

where x refers to the preferred direction and A is for the tensor polarizations, -|-, x. Note that a • denotes the inner 
product between a tensor mode, associated with the polarization of the gravitational wave which has two indices, 
and two vectors. For example x • • k = Xi kj . We should emphasis here that the above shape is the leading shape 

for the TSS. There are also some other terms which contribute in TSS. However, as we will show, their contribution 
would be sub-leading. 

As we will discuss, in models of anisotropic inflation there will be quadrupole anisotropy in two point function 
induced from the preferred direction so this TSS would be added on the top of that anisotropy. We try to use this 
anisotropy as a way to probe the tensor mode as well as to find the preferred direction. This can be done by writing 
the optimal estimator for GW. As we will show, there is a window in our parameter space which can yield signals 
detectable by Euclid and/or 21 cm surveys. 

The rest of this paper is organized as following. In Section |IT| we present the setup of anisotropic inflation in some 
details. In Section jlH] we calculate the TSS correlation function. Next, in Section EYi we consider the observational 
consequences of these analysis for the galaxy surveys. We end with Conclusion and leave some technical details into 
appendices. 


II. ANISOTROPIC INFLATION 

In this section we briefly review the setup of anisotropic inflation. For related works on anisotropic inflation see 
[30H121 , for a review see [331133] ■ 


A. The background 


Our anisotropic inflation model contains a 17(1) gauge field with non-zero amplitude at the background level. 
However, as it is well-known, the Maxwell theory suffers from the conformal invariance in an expanding backgrounds. 
As a result, the background gauge field A^it) is diluted exponentially during inflation. Moreover, the quantum 
excitations of the gauge field 5A^(t,x) during inflation will not be scale-invariant. Therefore, it is essential to break 
the conformal invariance. The easiest method to break the conformal invariance is to couple the gauge field to inflaton 
field to develop a time-dependent gauge kinetic coupling. With this discussion in mind the action of our system is 


S= 


- U(^) 


( 2 ) 
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where Mp is the reduced Planck mass while f{4>) represents the time-dependent gauge kinetic coupling which couples 
the inflaton field (j) to the gauge field. As usual, is the gauge field strength which is given by 


— dpi An — dnAfi. (3) 

The model of anisotropic inflation based on the above action was studied at the background level in and its 
analysis for curvature perturbations were performed in |31[ I34H40] . 

Note that in Maxwell theory /(^) = 1. However, as mentioned above, in order to break the conformal invariance, 
we need a time-dependent gauge kinetic coupling. As it has been shown in [501 - 140) . one requires / oc a~^ in order 
to obtain a scale invariant power spectrum for the gauge field quantum fluctuations. With this choice of the gauge 
kinetic coupling one obtains a constant background electric field energy density during the inflation. 

We choose the convention that the background gauge field is turned on along the x-direction so = (0, A 3 ;(t), 0,0). 
With this choice for the gauge field, the background space-time still poses the symmetry iny — z plane. The background 
is in the form of Bianchi I universe with the metric given by 

ds^ = —dt^ + a(t)^dx^ + b{tY (dy^ + dz^) 

= -dt^ + . (4) 


Here d measures the average of the Hubble expansion rate while measures the level of anisotropic expansion. 
However, the observational constraints on the level of anisotropy in curvature perturbations is not more than few 
percent [35] so ct/q; 1. As a result, one can treat the analysis perturbatively in this model. 

The background fields equations have been studied in [30|. It is shown that with the general gauge kinetic coupling 


f{(j)) = exp 




(5) 


the system reaches an attractor regime for c > 1 in which the ratio of the gauge field energy density pA (i.e. the 
electric field energy density) becomes a small but constant fraction of the total energy density p. It is convenient to 
define the fraction of the gauge field energy density to total energy density via 



2V 


As demonstrated in m, during the attractor phase one obtains 


R = 


I^H 

2 


( 6 ) 

(7) 


in which I = is a measure of anisotropy in the system and ep = —H /is the usual slow-roll parameter. The 
bound on R is determined by the quadrupole anisotropy in curvature power spectrum as we discuss below. 


B. Perturbations 

In the following, we review perturbations in anisotropic inflation setup. We divide the analysis into three parts: the 
quadratic action, the second order interaction terms as well as the third order interactions. 

In general, the perturbations in anisotropic setup originate both from the metric sector as well as from the matter 
sector. Specifically, the metric sector has the non-dynamical degrees of freedom dgo/x which should be integrated out in 
order to calculate the dynamical action. After integrating out these non-dynamical metric degrees of freedom one 
obtains many new terms in the dynamical action. However, as it has been verified in |34j . the perturbations from 
the metric are either slow-roll suppressed or would cancel with each other. Therefore, to leading order in slow-roll 
parameters and the anisotropic parameter R, the leading terms in the interactions come entirely from the matter sector 
perturbations. In other word, we do not have to consider the perturbations from the non-dynamical degrees of freedom 
in the metric which simplifies the analysis considerably. Technically speaking this corresponds to the decoupling limit 
in which one can neglect the gravitational back-reactions with the expense of inducing errors at the order of slow-roll 
parameters. 

To simplify the analysis further, we choose the flat gauge in which the curvature perturbations is given by the 
inflaton perturbations ^ = —^Scj). As a result, the metric takes the following simple form. 


= a{r])^ {—dif + -I- hij] dx^dx^) 


(8) 
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where t] represents the conformal time drj = dt/a{t). As mentioned above, since the leading corrections come from the 
matter fluctuations, one can safely set a = b] the errors induced with this assumption are higher orders in R and en- 
The perturbations hij above represent the tensor modes which are transverse and traceless: dih^j = 0 and ha = 0 
where the repeated indices are summed over. We denote the two different polarizations of the metric by /ix and /i+. 
Finally, for the gauge field perturbations, we choose the Coulomb-radiation gauge where 6Aq = 0 and diSAi = 0. 

We quantize the curvature perturbation, the gauge field perturbations and the tensor perturbations as usual. Let us 
start with the curvature perturbation 


^ Mpy/2eH 


(9) 


Going to the Fourier space, we can expand the curvature perturbation in terms of the annihilation and the creation 
operators a(k) and a'f(k) as follows 

/ d^k . ^ ^ 

(^7r)3/2 e'*'''‘C(k,??) , C(k,??) = C(fc,’7)a(k)+ (10) 

where, as usual, the creation and the annihilation operators satisfy the commutation relation [a(k), al(k')] = (k—k'). 

In the limit that we neglect the gravitational back-reactions, the wave function of the curvature perturbation 
corresponds to the profile of a massless scalar field in a dS background, 


Ck{v) 


Mpy/2eHk V k-q) 


( 11 ) 


Here we have assumed the Bunch-Davies vacuum. It is also possible to relax this assumption, as it has been considered 
in [45l|46]. However, in the following, we avoid this complication and simply assume that all fields start with the 
Bunch-Davies vacuum deep inside the horizon. 

The curvature perturbation power spectrum is given by, 

(C(ki)C(k2)) = (2^)35(3)(ki-fk2)Pc(A:i) , Vc^^Pdki) (12) 


The isotropic power spectrum, in the absence of anisotropy, is 


p(0) ^ 


8tt‘^€hMp 


(13) 


We calculate the anisotropic power spectrum induced from the gauge field later on. 

To study the gauge field and tensor perturbations we have to specify the wave number of the modes as well as the 
polarizations of the vector and the tensor fields. Although at the two point level it is enough (as usually assumed) to 
take the wave number vector in the x — y plane but this is not the most general choice at the third order level. The 
most general wavenumber is given by 


k = k cos d sin , sin 0 sin , cos (fj , (14) 

where (j) is the angle of the wavenumber with respect to the z direction and 9 represents its angle with respect to 
the X direction (note that the ^-direction is the preferred direction in our setup). Now, with the above form of the 
wavenumber, the vector perturbations are represented by 

MKv) = ^ , A^\)iKv) = A(^\)ik,v)a(\dk) + Addk,q)a4^i-'k), (15) 

A^l,2 

in which a(A)(k) and a|^j(—k) are the annihilation and the creation operators associated with the gauge field 

perturbations satisfying the commutation relation [a(;,)(k), a|j^,j(k')] = (5»,/ (5(^)(k— k'). In addition, e-^^(k) represents 

the gauge field transverse polarization satisfying fcie-'^^(k) = 0. With the above choice for the wavenumber, the gauge 
field polarization is given by 

e = (sin d, — cos 0,0) , e = (cos 9 cos sin 0 cos ip, — sin (/?) 


(16) 
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In addition we define A(i) = A and A( 2 ) = B. 

Similar to scalar field perturbation, in the absence of gravitational back-reactions, the quantum fluctuations of the 
gauge field perturbations are given by 


Akijl) = Bkiv) = 




1 - — 1 
krj^ 


(17) 


Now we consider the tensor perturbations. We decompose the tensor perturbations in terms of two polarization 
bases e(k) with s = x, -|- satisfying 


If (k)e ■ ^ (k) = Sss' , eff (k) = (-k), 


(18) 


in which a * represents the complex-conjugation. In addition, the traceless and transverse conditions ha = hijj = 0, 
yields 


eii^(k) = 0 , %eg^(k)=0. 


(19) 


The quantum operators hij (k, rf) in Fourier space can be represented in terms of their annihilation and creation 
operators by 


hy(k,r7)= ^ h(s)(k, 77 )e^®'(k) , hs(k, ry) =/i(s)(fc, ?7)as(k)-b/i^,)(fc, r7)4(-k), 

S = + , X 

with the commutation relations [as(k),aj(k')] = (5ss/5(^)(k— k'). 


( 20 ) 


We also need to fix the tensor polarization. With the wavenumber given by Eq. (14), the most general tensor 


polarization is given in Appen. However, it turns out that for the calculation of the TSS correlation we still have a 
freedom to choose one of the fields in the x — y plane. So in order to simplify the analysis, we can simply assume that 
the tensor mode is completely inside the x — y plane corresponding to choosing ^ in the following. With this 
choice the polarizations efj(k) and efjCk) become 


^ sin^ 9 — sin 9 cos 9 0 

efjCk) ( — sin0cos0 cos^ 0 0 


1 


0 


0 


-1 


1 


0 0 — sin 9 

ef (k) I 0 0 cos6» 


- sin 9 cos 9 0 

As usual, the profile of the tensor excitations has the standard massless form 

Correspondingly, the power spectrum of the tensor perturbations is given by 

(77ki)77k2)) = ( 2 ^) 373 ) (ki + k2)P7fci) 




In the isotropic case, the tensor power spectrum has the standard form 


p(0) ^ 


2H^ 

7r2M|, 




(0) 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


Therefore, defining the tensor-to scalar ratio r = we obtain r = 16e// for the isotropic theory. 

C. The anisotropic power spectrum 

As discussed before, the presence of a massless vector field interacting with the infiaton field induces anisotropy in 
scalar and tensor power spectrum which can be detected observationally. All interactions responsible for generating 
anisotropic power spectra are pre-dominantly generated from the matter sector Lagrangian [34] 




(26) 
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The anisotropic power spectrum is usually parameterized via m 

n (l+5*cos2 0fe) , (27) 

in which 0^ is the angle between the preferred direction n (in our example n = x) and the wavenumber k and g* is the 
amplitude of quadrupole anisotropy. Using either the standard in-in formalism, as performed in [SI1IM1I3H1, or the 
alternative 6N formalism |d5l Kfti] , the amplitude of quadrupole anisotropy in the attractor limit of m is calculated to 
be 


g, = - 24. (28) 

Using the Planck data, the bound |g*| < 10“^ is obtained in [32]- As a result, with Ne ~ 60 to solve the horizon 
and the flatness problem, one concludes that / < 10“^. Translating this to parameter c appearing in gauge kinetic 
coupling, Eq. ([^, one concludes that c — 1 < 10“^. On the other hand, in order for the system to reach the attractor 
regime during which the gauge field energy density is given by Eq. Q, it is assumed in [31] that (c — l)iVe > 1. We 
clearly see that this condition is not consistent with the observational bound on g* [45]. It is therefore necessary to 
revisit the anisotropy power spectrum without relying on the assumption of attractor regime as performed in |48j . 
In the limit where one neglects the slow-roll corrections and assuming the total number of e-folds N^. are finite, the 
appropriate limit of the results in [45j yields 


9 * = -24 


2R 

(-H 


Nf 


(29) 


We have specifically checked that this result is also obtained using the SN formalism [49]. We note that Eq. (29) is 
consistent with the attractor limit in which from Eq. 0 one has 2R/€h = I- 

In conclusion, as we have checked in |49j . the previous results of anisotropic curvature perturbation power spectrum 
and bispectrum obtained in [811 l84H4()j go through directly without relying on the attractor assumption upon replacing 
/ —>■ 2R/eH- In particular, the amplitude of quadrupole anisotropy scales as which is supported from the 
accumulative contributions of the super-horizon IR modes which have left the horizon and becomes classical |38j . 


III. THE TSS BISPECTRUM CORRELATIONS 


Having presented the free held mode functions, now we are in the position to consider the interaction Lagrangians 
and calculate the TSS bispectrum correlations. As mentioned before, the leading interactions for anisotropies and the 
TSS correlations comes from the matter sector Eq. (26). 


Expanding the above action around the background values, one can obtain the interaction Lagrangian at the second 
and the third order levels. As usual, one can start by writing the interaction Lagrangian in the real space, as presented 
in Appendix]^ and then try to transform it into the Eourier space. It turns out that the Hnal result for our most 
general choice of wavenumber k is very complicated. While we present the full results in Appendix [^ in what follows 
we only consider two simpler limits of the system. As we will see, these simple examples will be useful for us to better 
understand the leading order corrections. 


A. The all planar shapes for the bispectrum 


Here we consider the simple case in which all wavenumbers are co-planar extended in x — y plane, corresponding to 
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1. The second order interaction: 


In this case, the second order exchange vertex Lagrangians in the Fourier space have the following form, 


= -QV2RMp j d^ksiT? 9ka^ {-vi) ^ Ck/»+k 

= —‘iMpV&R J (f k sin 9k — ^ 


Lh+A = MpV^R d^ksin9k ) A{^h. 




k^+k 


Lhy,B = MpV^R J d?ksin9k ) Bk (/ixk) 


fa 


(30) 

(31) 

(32) 

(33) 


where 9k is the angle between the preferred direction n and the wavenumber k (in our example h = x) and so on. 


2. The Third order interaction: 


Here we present the third order action. Since we would like to calculate the TSS bispectrum, we only present the 
leading contributions into the TSS correlation. As we will see, at the leading order, there are five different diagrams 
which contribute into TSS correlation. They correspond to the following third order terms in the cubic action, , 
rSSV^ l^SVT^ l^svv where S', V and T respectively refer to C, A and hp. 

It is very important to note that beyond the above terms there are some other terms at the cubic action which turn 
out to be smaller than the leading corrections that we are interested in so throughout this work we do not consider 
them. Fortunately there is an intuitive factor in order to determine the leading terms which we specify in what follows. 
Moreover, there are also many other terms inside the above five terms which either do not appear in TSS or lead 
to a sub-leading contribution as compared to the leading term. We also neglect these terms and only mention the 
potentially leading contributions. 

Now, before going ahead, it is worth to introduce the other terms inside the above five terms and mention briefly 
why we are not interested in them. There are two types of terms in this regard: 

• Terms which contain the spatial derivatives of the gauge field corresponding to the magnetic field. As has been 
shown in the Appendix |D 2[ the final result for this part is logarithmically suppressed compared to the electric 
field coming from the time derivative of the gauge field, so these terms do not contribute at the leading order. 
Intuitively, this comes from the fact that magnetic field decays outside the horizon while electric field is frozen 
on super-horizon scales. 

• Terms which either contain B or h-^. However, in the planar case, y: = f, C only interacts with the scalar fields 
such as A and hp at the second order level. As a result there is no leading sources for these types of interactions. 


Thanks to these simplifications for the planar case, the third order interactions are simplified substantially given by 

^2 ' 


L 


SST 


-SSV 


-12V2RMp j dfp J d^qsin^ 6i(_p_g)CpCq^-H(-p-q) 
-SVfkRMp J d^p J d^q sin 9qC(^_p_q^CpAq 


L^VT ^ 2Vl2i?Mp J J cl^'fsinSpCOS (0p - 0q)C(_p-q)/i+pAq 

rSW ^ y2 J J d39COs(0p-0,)C(-p-q)A'pA'q 

^VVT ^ _y^f2 J J d^q cos {9p-9(^ppq)) cos {9q-9(^ppq))hp^_p_q)A'pA'^ 


(34) 

(35) 

(36) 

(37) 

(38) 






(A) 





(D) 


(E) 


FIG. 1: Feynman diagrams for the Tensor-Scalar-Scalar cross-correlation with different types of exchange vertices. 


B. TSS correlation function for the planar shape 


Having presented the cubic interactions, in what follows, by using the in-in formalism, we calculate the TSS 
bispectrum for each of the five mentioned pieces in the planar limit. As we shall see, this example turns out to be 
extremely useful for us in order to determine the leading term in the full TSS bispectrum. 

However, before going ahead, we emphasis that the standard in-in formalism is written in terms of interaction 
Hamiltonian. So in order to use it we need to calculate the Hamiltonian from the interaction Lagrangian of the system. 
Moreover, in the presence of the kinetic interactions, Hint is not simply —Lint- However, as it has been shown in EZ], 
the leading corrections in the correlation functions come entirely from the substitution Hint = —Lint and we can safely 
use this prescription in the following analysis. 

In order to have a better intuition about the potentially important terms it is worth to take a look at the Feynman 
diagrams in Fig. representing how the above five terms can contribute into the TSS bispectrum. From the plot we 
see that, and correspond respectively to the cases A, B, C, D and E in Fig. 

Now, before doing any analysis, an important question is how can we determine the potentially important contribution 
in the TSS bispectrum? In order to find the answer, let us take a loo k at the simplest possible contribution in TSS 
which comes from represented by diagram A in Fig. From Eq. (341 we see that this interaction is proportional 

to i? ^ 1 which can be used to determine the potentially important terms in the TSS bispectrum. As we will show, all 
terms in Fig. are proportional to R. 

The next question, related to the above question, is whether all diagrams presented in Fig. [l] are important in 
performing the in-in analysis. Actually, it turns out that not all of these terms are the same order. After performing 
the time integrals the leading contributions come from corresponding to diagram E in Fig. Having this said, 

in the following, we calculate TSS bispectrum for all diagrams in Fig. Here we only present the final result for each 
diagram and leave the detail analysis to Appendix. 


• SST: This corresponds to diagram A in Fig. [^for which we need to calculate in-in integral to first order. The 
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result is obtained to be 

(^+kiCk2Ck3 iVe)) 




f kf + k2 + 
[ kl 


-P/i+ki^Ck2 + (1^2 ^ ka) , 


(39) 


where we have defined the number of e-folds as = — In and is the angle between the preferred 
direction h and the wavenumber ki (in our example h = x). 

• SSV: This term corresponds to diagram B in Fig. In this case, we have to calculate the in-in integral to 
second order. The final result is obtained to be 


(^-('kiCk2Ck3 iVe)) 



{mi - 2Ne) sin^ Ok, 




^);k2^/i+ki + (k2 ^ ka) 


(40) 


• SVT; This term corresponds to diagram C in Fig. In this case, we have to calculate the in-in integral to 
second order, yielding 


(h+k.Ck2Ck3(?7e)) = {m!-2N,) 


kf + fcf 


sin6»fcj sindfc^ cos (0^^ - j -Pckak^ + (k 2 ^ kg) 


• VVT: This term corresponds to diagram D in Fig. For this case, one has to expand the in-in integral up to 
third order in the perturbations yielding 


(^-rkiCk2Ck3 (^e)) — 


(“S?) 


Nt-Ni + iNe) (§ 

^3 


sin 6 »fc 2 sin 6 »fc 3 cos {Ok^ - cos (dfcj - dfc 3 )Pck 2 ^’?i+ki 


+ (k 2 ■!-> ka). (41) 


• SVV : This term corresponds to diagram E in Fig. Like in the previous case, one has to calculate the in-in 
integral to third order yielding 




\ en J 

sinOk^ sin 6 »fc 3 cos {Ok^ - 6 »fc 2 )Pc(k 2 )-P?i+(ki) + (k 2 O kg) 


(42) 


Now adding the above five contributions the final result is obtained to be 


(^+kiCk2Ck3 (?7e)) — (^+kiCk2Ck3 (?7e))y4 + (^+ki Ck 2 Ck3 (^e))^ + (^-l-ki Ck 2 Ck3 (^e))c 
+ (^+kiCk2Ck3 (?7e))£) + {k+WiCk^CksiVe))E 


\ 

( X • e+(ki) • X - (kg • x) (^x • e+(ki) • kg) 


sindfc^ sin 6 »fc 2 cos {dk^ - 0 k 2 )Pci'^ 2 )Ph+{'i^i) + (k 2 kg) 


PckisPh+k^ + kg o k 2 j , 


(43) 


where we have only considered the leading contribution in the final result which comes from SVV represented by 
diagram E in Fig. 
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C. Planar tensor and non-planar vector fields 

As a second example, in the following, we consider the case where the tensor held is in the x — y plane, i.e. i^tensor = f, 
while the polarization of vector held has the general non-planar form. In this case, the situation is more complicated 
than the hrst example studied above and the above 5 diagrams are converted into 13 diagrams. As a result the situation 
is much more complicated than the previous example in which all modes were assumed to be coplanar. However, one 
does not need to worry about these complexities. This is because, as we have already proved in the previous example, 
there is a hierarchy between different terms and thanks to this hierarchy, we only need to consider the last diagram, 
the diagram E in Fig. Since this diagram comes from the SVV vertex, in the following, we only need to consider 
this vertex and do not consider the other 12 diagrams. Having this said, we only consider the second and third order 
interactions relevant for this diagram. 


1. The Second order interaction: 


The relevant second ordered interactions are given by. 


Lsv 

Lvt 


—4:MpV2R J ^ (sin0fcAj^-I-cos0fc cosyifcH^) 
Mp\fiR J ^ (sindfeA'^k^+k + sindfeH^i^/ixk) ■ 


(44) 

(45) 


2. The Third order interaction: 


The third order interactions is given by. 


Lsvv = 2/^ J d^pd^qC(-p-q) 


cos (Op - d,) ApAq -I- 2 sin (dp - d,) cos (^pAqHp -I- sin i^p sin ipqBpB'^ 


(46) 


D. TSS correlation for the planar tensor and general non-planar vector fields 


Here we calculate the TSS correlation function for the planar tensor field with the vector field having the general 
non-planar wavenumber. As we mentioned above, we consider only the most relevant diagram for this case which comes 
from the SVV type interaction. As it has been shown in Fig. in this case there are three different contributions that 
are worth to be considered. 

♦ First we calculate (h+kiCk^CHiVe)) which is 


(^+kiCk2Ck3 (de)) — 


( mV2R' 


n , Vp 


y £h J \ '^ “ 3 

sindfei sindfcj cos (d^^ - d^j) -hsindfe^ cosdfcj cos^ sin(dfc3 - d^J) Pcks^/i+ki + (ks ^ ka) 


- 25*^6 


(x ■ e+(ki) • X - (ks • x)(x • e+(ki) • ks)^ PcksPh+k^ + (ks ^ ka) 


(47) 


where, as before, x refers to the unit vector in the preferred direction and we have used Eq. (29) to express the final 
result in terms of the amplitude of quadrupole anisotropy (/*. 


♦ Now we calculate (dxfciCfc 2 Cfc 3 (de))i obtaining 


(hxkiCk2Ck3 (de)) — 


f my/2R' 


1 


n: - n: + -N, 


sindfej cosdfcj sinvJfca cos + (ks ka) 


y £h J \ '^ “ 3 

(^x • e'^(ki) • X - (ks • x)(x • e^(ki) • ks)) Pc('k3)Ph^ (ki) -h (ks O ka) 


3 


(48) 
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FIG. 2: Leading order corrections on the TSS cross-correlation function for the non-planar case . 

There are two non-trivial features in these TSS correlation functions that in the following we elaborate in some 
details. 

• First of all, we emphasis that there are two different sources for the TSS correlations. One comes from the isotropic 

part of the action, which is similar to the single field inflation and is almost canceled by the other terms. The other 
term which we are interested in here is indeed due to the anisotropic matter corrections to the FRW universe which is 
determined by the amplitude of the quadrupole anisotropy in power spectrum 5 * = — One feature of this 

correction is that it has primordial IR origin and it can not be canceled by any late time effects such as the non-linear 
effect as well as the projection effect. Yet, this correction is under control as it is proportional to the quadrupole 
anisotropy at the two point level. 

• The second interesting point is that the shape of the TSS is not the same as in the well known fossil shape given 
in [ 2 ]. This is because in our model, in addition to the wavenumbers and the polarization tensor, we also have the 
effect of preferred direction, i.e. the statistical anisotropy. As a result one expects that this preferred direction would 
show up in the final results. 


IV. OBSERVATIONAL SIGNATURES 

In this Section, we present different observational signatures of our anisotropic set up. 


A. SS vs TSS power asymmetry 


As we have seen breaking the statistical isotropy leads to a quadrupole anisotropy at the two point level of the 
primordial curvature perturbation with the following form, 

Pc(Ks)=Pc°(i+5*(Ks-S)') , (49) 

with Qt = — ^e- Here we have shown that in addition to the above quadrupole anisotropy in the power 

spectrum, there is also another prediction from anisotropic inflation which comes from the anisotropic correlation 
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between the short scalar modes and long tensor modes, 


^(Ks) 


'«(A)(Kl) 


Pciks)A^^h^\Kj^), 


(50) 


where is given by, 




^144i?A3 j 
= 35*A^e( x'x-’ 


(^x*xJ-(Ks-x)(x'Ki)) 

-(Ks-x)(x'Ki)) 


(51) 


So before going ahead, it is suitable to compare the results from the two-point level with the result at the cubic order. 
We can compare the amplitude as well as the shape of these two functions with each other. Let us first compare the 
amplitudes. The ratio between the amplitudes of the power asymmetry at the two point as well as the three point 
levels is given by. 


<5(Pc(Ks)U(„(ko) 

<5(Pc(Ks)) 


= 3Aeh(;,)(KL). 


(52) 


This means 


For reasonable values of total number of e-foldings, say Ne < 100, the above ratio would be of order 10“^. 
that the quadratic corrections are more important than the cubic corrections in this set up. 

As for the shapes, comparing Eq. (49) with Eq. (50), we see that the shapes of these asymmetries are not the same. 


So in principle one can distinguish them from each other. 


B. Quadrupole Anisotropy in the Galaxy Power Spectrum 

Here we hnd the shape of the galaxy power asymmetry within our anisotropic setup. In order to do this, we would 
use the orthogonality of the spherical harmonics, Tlm- In addition, it is clear that the asymmetry is in the form of 
quadrupole which means that we are only interested in the L = 2 case. So we would have, 

/d^fcPc(|fc|)A^^hg)(KL)T;M 

Jd^kPci\k\)Y*, 

J d^kP^(\k\)A^^h[^{K^)Y*j^ (53) 

We note that the quadrupole components are not a good observable by themselves and in order to get the observables 
one needs to calculate the root-mean-square of the above quantity which is given by, 

= ( E \Q^mA- (54) 

' M=-2 ' 

In addition, as we shall see below, depending on our purpose for the detection of the tensor mode, we may also need 
averaging over the long mode and in the case that we do not have any information about the preferred direction, we 
also perform an averaging over this direction. Keeping these in mind, in the remaining part of this section, we try to 
calculate the root-mean-square of the quadrupole moment. As soon as we calculate it, depending on our purpose, we 
may also do the averaging over the long-mode as well as the preferred direction. 

After an extensive analysis, which is presented in the App. (H, the root-mean-square of the quadrupole is obtained 
to be 


Q 2 M = 


q" = (Stt) {g^N^f Pu{Kl) (l - (Kl • x)^) A - (dOTr^) ^ ^ F* (x)r*, (x)Ff,(KL)yiV (Kl) 


2 1 1 
000 


2 1 1 


2 1 1 


t q s J \ —t q s 


(55) 
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where 


I I, U 


is the Wigner 3-j symbol. 


m mi m 2 

Averaging over the long mode ; 


In the next step we calculate the average of the root-mean-square of the quadrupole anisotropy over the long mode 
as follows 


Q \l ^ 


rKs„ 




d^KLQ (Kl,x) 


= {8Tr){g^N^f 


H V Ksrmn \ 


^Mp J / 

s,s' ,q,q' — —l t— — 2 


(i-(Kl-x)2) (^ 1 - (dOTT^) 


2 1 1 
0 0 0 


= (Stt) {g^NeY ( ^ 


In 


f Ksr, 


2 1 1 
0 0 0 

2 1 1 
0 0 0 


2 1 1 
t q s 

2 1 1 
t q s 


\ ^Lmin 

2 1 1 
—t —q s' 

2 1 1 
-t q' s' 


Stt f SOtt^ 

"3 


+1 


2 1 1 
t q s 

+2 


2 1 1 
-t q' s' 


V 3 


s,s' ,q—— 1 t— — 2 

)\/J E E E (-ir+’Vi*(x)yi*,(x)y2-;„(x) 

s,sYq,q' — — lt— — 2 m——2 


\ 3 

/647r3\ /45 


2 1 1 
m —q' —q 


(56) 


Interestingly we see that in this model the averaged power quadrupole anisotropy is logarithmically sensitive to the 
number of inflationary modes in the survey. 

Averaging over the preferred direction: 


So far we had assumed that we know the preferred direction in the sky (the direction of anisotropy) from other 
experiments. Now one can ask what will be our observables when there is no knowledge about the direction of 
anisotropy. In order to calculate the observed quadrupole anisotropy in this case one needs to average over the preferred 
directions. So in the following we also calculate the averaged power quadrupole over different preferred directions. We 
start with Eq. (56) and perform averaging over different values of the preferred direction 


Q^\f = J 


= (Stt) (g^Ne) [ jp; j 


dd Y , ( Ksmin \ ( 327r2 


2 1 1 
0 0 0 

2 1 1 
0 0 0 


2 / \ 2 
^2 1 1 ' 


t q s 

2 1 1 
t q s 


\Kl min ) 

-h (240-2 


+1 


SOtt^ 

+1 


+1 +2 

E E (-!)«+■ 

s,q——l t— — 2 


E E (-1)’*'" 

s,s',g',g' = —1 t,m— — 2 


2 1 1 

-t q' s' 


2 1 1 
m —q' —q 


2 1 1 

—m —s' —s 


(57) 


C. Optimal Estimator for the Anisotropic Tensor mode 

Since the TSS correlation function is proportional to the tensor amplitude, one can use this set up to estimate the 
tensor modes within the anisotropic inflationary models. In the following, we try to do this in our set up. We start by 
writing down the minimum variance estimator for the tensor mode and then will continue with calculating the optimal 
estimator. Then we discuss briefly about the observability of the tensor modes in this set up by using the current and 
future galaxy surveys. 

As it has been mentioned in [5], in the presence of the tensor-scalar-scalar correlation, each pair of (5(k2) and ^(ks) 
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provide an estimator for the tensor mode with K = k 2 + ka, where we have redefined (— ki) as K. 


hix){K) 


6{k2)S{ks) 

P{k2) 


(A*^(x,k2)e(f(K)) 


(58) 


In the following, we assume the null hypothesis, in which in the absence of the tensor mode the density perturbations 
are statistically isotropic and Gaussian. Of course in our current set up it is not completely true because as we have 
mentioned above we have broken the statistical isotropy so we have some anisotropic corrections in the power spectrum 
of the density perturbations. However, the level of anisotropy is very small and one can neglect them at this level. 
Naively, it means that we can write ((5(k)(5(k')) = where V appears here due to the fact that we have 

discretized the density function and is the Kronecker delta function and = P(k) + p("^(k) is the measured 

matter power spectrum, which is including both the signal, H(k), and the noise, p(")(k). 

We can then calculate the variance of the above estimator, 


Variance = 2VP‘°*(k2)P‘°‘(k3) 




(x,k2)el^)(K) 



(59) 


where the factor 2 comes from the permutation and k 2 ) = (x, k 2 )P(k 2 ). 

The minimum-varince estimator for the h(^x) (K) can then be calculated by summing over the whole individual pairs 
and with the inverse-variance weighting, 


V)(K)=pj:)(K)^ 

k 


Sb(x,k,k-k)eif(K) 

2yptot(k)ptot(K-k) 


<5(k2)<5(k3) , 


(60) 


where in order to emphasis the symmetric beha vior under k^ O k 3 , we have presented this function more precisely as 
S®-' (x, k, K — k). The reason is that, from Eqs.(47) and (48), it is clear that due to the permutation between k 2 and 
k 3 the answer should be symmetric. So far we did not take this into account because we are in the squeezed limit 
and these two modes are equal to each other at the leading order. However, it is important to keep in mind that the 
precise notation is the one given by Eq. (60). 

In addition the noise power spectrum is also given by 


4”) (K) = 


E 

L k 


'(x,k,k-k)e(;)(K) 


1 -1 


2V- k\) 


(61) 


In the most general case, the tensor mode, h^x) (K) arises as the realizations of the random fields with the power 
spectrum, P(;,)(K) = A(x)P^x)(^)^ where A denotes the polarizations of the tensor mode while refers to the 
amplitude of the tensor mode and P^'^^(K) is the fiducial power spectrum. 

Now each Eourier mode of the /i(;),)(K) provides an estimator for the amplitude of the tensor field as follows 


Aa)(K) = 


^(a)(K) 


1 -1 


V- 


h(x)(K) 


-^(a)(K) 


(62) 


where we have subtracted the noise contribution to unbias the estimator. Now, as we mentioned above, each pair of 
the density fields would act as an estimator for the gravitational wave. So if we have a huge number of the density 
fields then according to the central limit theorem h(x) (K) would be very close to a Gaussian variable. As a result the 
variance of the estimator in Eq. (62), under the null hypothesis, would be 


Variance of Amplitude = 2 


^/a)(K) 


-2 




(63) 


Now adding the estimators for every Eourier modes with inverse variance weighting yields the optimal estimator^ 


P^) = alY, 


^/a)(K) 
^(a)(K) 


K.A 2 


V- 


h(x)(K 


-^^(a)(K) , 


(64) 
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where we have 


= E 


P(x)i^) 


if,A 2 


4")(K) 


(65) 


Comparing the above optimal estimator, Eq. (641, with the minimum estimator, Eq. (60), we see that the correlation 


of the density fields with the fossil tensor modes appears as a four point function in terms of the density helds. Now, 
according to the above estimators, a signal is detected if the null hypothesis estimators defined above is found to 
depart at > Scr from the null hypothesis. 

Now we evaluate the smallest amplitude Ah which can be detected within a given survey. In order to do this, we 
assume that the fiducial power spectrum is nearly scale-invariant, say C/t(K) = AhK^''~^ with Uh — 0. We also use 
our results for the squeezed limit of the fossil-density-density field bispectrum, Eq. (501. We then take the continuous 
limit of the summation in Eq. (61), i.e. by using —>• V f d^k/{2Tr)^. We also assume that for k < kmax we have 

{-jMht) ~ 1, where k^ax is the largest wave number for which the measured power spectrum would have a large signal 
to noise, and for any k > kmax this ratio is equal to zero, i.e. we are dominated by the noise. Finally, we also have to 
take into account the direction-dependence of the tensor field. We consider the most generic case, given in Appendix 
A, and investigate the observability of a signal for this generic case. Then at the end of the day in order to calculate a 
we need to perform averaging over different directions of GW. With this discussion in mind would be. 


(pf)(K)) 


+ 

( — 28 cos'^ 9 COS^ (j) + 29 COS^ 9 COS^ (j) ~ COS^ 0 -I-1) . 


-1 ^ ( gjN? 

\ dOTT^ 

(glNf 


km^^ (7cos^^ ^cos"^ (j) + Idcos"^ 0cos^ (j) — 13cos^ 9cos^ (/) -I- 7cos^ 6 — 15cos^ 0-1-8) (66) 


(67) 


We can then plug back the final expression for the noise power-spectrum into Eq. (65) in order for the tensor amplitude 
to be detectable at the > 3a, 


= 



( 68 ) 


So in the presence of the tensor mode with the amplitude = 3 x 10 (corresponding to tensor to scalar power 
spectra ratio ~ 0.1 ) in the matter two point function, kmax/kmin would be 


kmax/kmin = 12800 (g,W)^ ■ 


(69) 


For 5 *^ 2 and ^ 8 we obtain respectively kmax/kmin ~ 8000 and kmax/kmin ^ 3876. These predictions seem 

to be detectable in the galaxy surveys like EUCLID and/or the 21cm surveys. 


D. Estimator for the Preferred Direction 


So far we had assumed that we know the direction of anisotropy in the sky from another experiment. However, 
it might be that there is no knowledge about the direction of anisotropy in the sky. So in this case it is worth to 
determine the shape of the estimator that we tried to calculate. For this purpose, we start from Eq. (60) and construct 
the following estimator. 


h|™)(K) = J dx^h^)Y:^{^) 


• , _ A*Hx,k,K-k)e))'VK) 

dx2P("^) (K, x) ^ ^ (*) 


2V 


(70) 


Since in this case the direction of anisotropy is not known we have written explicitly the dependency of p|"^(K,x) on 
the preferred direction, x. 
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In the following, we try to approximate the above integral and hgure out which values of I would be non-zero. For 
this purpose, it is convenient to neglect the dependency of on the preferred direction. So we are left with 

the following integral, 

I (71) 

where we have neglected the denominator as well as the summation over k and the tensor. It is straightforward to 
show that the above estimator looks like a quadrupole. In order to show this, one should rewrite x* in terms of IdmS. 
Then according to the additional momentum theorem, we see that the hnal value of I can be only 0, 1 or 2. However, 

since the integral of the three spherical harmonics is proportional to ^ ^ ^ g ^, we see that for I = 1 this Wigner 

three-j symbol is zero. So we are left with either Z = 0 or Z = 2. Among these two different cases, the first one is 
completely constant so we can imagine that it contributes to trace of the e^- which is zero. So at the end we are left 
only with I = 2 which proves that the result indeed looks like a quadrupole, 

« ( 0 J J ) . (72) 

V. CONCLUSION AND DISCUSSION 

While the searches for the primordial GW in CMB are ongoing, it is worth to consider the large scale structure 
surveys for shedding light on the imprints of primordial GW. The goal of these surveys is to study the two point 
correlations of the density fluctuations. According to the null hypothesis, there is no off-diagonal correlation between 
different Fourier modes of the density fluctuations. However, in the presence of a long wavelength tensor mode different 
Fourier modes couple to each other. As a result this off-diagonal correlation function can be thought as a way to probe 
the primordial GW. On the other hand, since the key point in these surveys is the power spectrum of the density 
perturbations, there should be a hierarchy between the scale of the tensor mode, assumed to have a long wavelength, 
as compared with the typical scale of the density perturbations, supposed to have small wavelength. In another word 
one should consider the squeezed limit of the Tensor-Scalar-Scalar bispectrum. However, the tensor mode wavelength 
can be larger, comparable or smaller than the Hubble scale at any time. It has been proven that while in the first 
two cases the amplitude of GW is still nearly constant, it is not the case for the last case, i.e. when the tensor mode 
re-enters the horizon. In this case, it starts oscillating and decay but even in this case its effect survives even after the 
decay of the tensor mode [531. That is why this effect is called the fossil effect. Moreover, the shape of this off-diagonal 
correlation function looks like a quadrupole anisotropy. 

There are three different contributions in the above TSS coming from the inflation, the non-linear mode coupling 
during the radiation and the matter domination era and finally the projection effect. It has been shown previously that 
for single field slow-roll inflationary models, there is a cancellation between the late time projection effect and both of 
the primordial inflationary contribution as well the non-linear mode coupling in the super horizon and sub horizon 
limit respectively. Therefore, the final correlation function is small in these models. However, the above cancellation 
does not exist for the models beyond the SFSR inflation and thus the quadrupole anisotropy can be large enough in 
these models which can be also thought as a way to distinguish between the single field and multiple held models. 

On the other hand, in models of anisotropic inhation in which there is a preferred direction in the early Universe and 
the three dimensional rotational symmetry has been reduced to a planar two dimensional symmetry, the primordial 
power spectrum has quadrupole anisotropy in the absence of any tensor mode. Therefore, it is worth to consider these 
models in detail and hgure out a way to distinguish between the predictions of these models and models in which the 
rotational symmetry is intact. 

In this paper, we have calculated the quadrupole correction in the power spectrum in models of anisotropic inhation 
constructed from the dynamics of a C/(I) gauge held and have shown that this quadrupole is completely aligned 
along the preferred direction. We then considered the three point TSS correlation function and determined the shape 
of the quadrupole anisotropy for this case as well. As we have shown, in this case the quadrupole anisotropy is 
partially aligned along the preferred direction. This means that the models of anisotropic inhation have a unique and 
distinguishable signature in large scale surveys compared to models of inhation based on scalar helds. 

Finally, we considered the observational consequences of this quadrupole anisotropy and used it as a tool to estimate 
the contribution of the primordial GW as well as to hnd the preferred direction. We have presented the optimal 
estimator for the amplitude of the gravitational waves and have shown explicitly that there is a window in the parameter 
space where there is a 3a detection limit on the amplitude of GW where the signal is detectable in the galaxy surveys 
like Euclid and 2I-cm. 
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Appendix A: The most generic tensor polarization 


In this appendix, we present the most general form of polarization for the tensor modes. Given a wave number as 


k = fc(cos0sin(^,sin0sini^,cosy)), the most general tensor polarizations eMk) and ef (k) are given by 


e+ 


e^- 
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1 — cos^ 6* — cos^ 0 cos^ — sin0cos0(l + cos^ (^) sin 93 cos cos 0 


— sin0cos0(l + cos^ (/?) cos^ 0 cos^ (^ + cos^ 0 — cos^ sin 0 sm (/? cos 

■ sin^ Lp 


sin ip cos p cos 0 

2 sin 0 cos 0 cos p 
— cos(/3(2cos^ 0—1) 
— sin p sin 0 


sm u sm p cos p 

— cos (/3(2 cos^ 0 — 1) — sin(/3sin0 
—2sin0cos0cos(p sin(/jcos0 
sin p cos 0 0 


(Al) 

(A2) 


However, as we mentioned in the text, still we have the freedom to assume the planar tensor modes for the whole of 
our discussions in this work, so we can simply put p = ^. With this choice the polarizations are given in Eqs. (21) 
and (221. 


Appendix B: Interaction Lagrangian in the real space 

In this Appendix, we present the third order interaction Lagrangian in the real space. Our starting point is the 
matter Lagrangian, 


= . (Bl) 

The goal is to calculate the third order interaction Lagrangian which can be obtained by expanding the above action 
up to the third order around a given background. Since the resulting Lagrangian contains a lot of terms, and some of 
them are irrelevant to our discussion here, we restrict our attention only on to the following five sets of interactions, i.e. 
p^SST^ ^SSV^ p^SVT^ j^SVV ^VVT^ addition, as we will show in the next Appendix, since the magnetic field is 
exponentially suppressed compared to the electric field, in the following we only consider the time derivatives and do 


not write down the spatial derivatives. 

lSSt ^ _2ARMI ax)ax)KAx) (B2) 

lSSV ^ -sVMMp C{x)C{x)SA[{x) (B3) 

jl^svT ^ q-ZeRMp C{x) {SA[{x)ha;x{x) + SA2{x)hxy{x) + dA'^{x)hxz{x)) (B4) 

lSvv ^ 2fax){{6A[{x)f + {6A'2{x)f + {6A',{x)f) (B5) 

L^VT ^ (^hxx{x){ 6 A\{x))‘^ + 2 hxy{x)SA[{x)SA 2 {x) + 2hxz{x)6A[{x)SA'^{x) + hyy{x){ 6 A 2 {x))‘^ 

+2hy^{x)6A2{x)SA'^{x) + h^^{x){5A'^{x)Y^ , (B6) 
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where we have used the following useful expressions, 


M; = MpVMi-r,)-h 


= Mp\l — 
sh 

if' 


w 

The remaining part is to go to the Fourier space which is standard and we do not present it here. 


(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 


Appendix C: Interaction Lagrangian in the Fourier space 

In the following, we present the second and third order actions in the Fourier space for the most general choice of 
the wavenumber as well as the vector and the tensor polarizations. 


1. Second ordered action 


Let us start with the second order action. 


Lsv = 
Lst = 
Lyr = 


—AMpVOR J ^ C-k (sin^fc^k + cos0fc cost^fcBk) 
—6V2RMp J sin^ 6lfe^_kh+k 
MpVZR J d^fc ^ ^sindfcA'_ij/i+k + sindfcB(^kdxk 


(Cl) 

(C2) 

(C3) 


2. Third ordered action 


Here we present the full third order interaction Lagrangian. 

.2 


L^ST ^ j d^p y'd3gsin2d,C(-p-q)Cp/i+q 

L^sv = -gy/^Mp J d^p J d^gC(-p-q)Cp ^sindgAq + cosd, cosi^gHq 
l^SVT _ 2\/l2RMp J J ^^dCi-p-q) ^sindpCOs(dp — dg)d+pAq 

+ sin dp cos(/?g sin (dp — 6g)hppB'^ + sin dp sint^qdxpBq^ 

j^svv ^ 2f J d^p J d^ qC(^^p_q^(^cos {Op-0q)ApA'^-2 sin {Op-0q) cos + 

X cos Pp cos pq + sin pp sin B'^B'^ 


cos (dp — 0q) X 


(C4) 

(C5) 

(C6) 

(C7) 
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L 


VVT 


J J ^(p+ 9 ) cos (^p + 6 >g) + sin 6 '(p+,) cos 0 (p+,) sin ( 0 p + 6 »,) + cos 0 p cos 0 ,^ 

x/i_|_(„p_q)ApAq + ^sin0p cos9q cos (fiq — cos^ ^(p+q) COS ifig sin {9p + 9q) + sin (20(p+g)) sin0p sin0q cos (pq 
x/i+(_p_q) ApBq — ^ cos^ ^(p+g) cos pp sin {9p + 9q) — cos 9p cos pp sin 9q — sin (20(p+g)) cos 9p cos ipp cos 9q 
x/l+(„p_q) AqSp + ^2 sinyjq COS {9p+q - 6»p)^ /lx(-p-q)^p5q + ^ COS 9p COS Pp cos 9q cos Pq - sin^ Pp sin^ Pq 
— cos^ ^(p+g) cospp COS pq COS (9p + 9q) — sin (20(p_|_q)) cos0p cos pp sin 9q cos pq^ /i+(_p_q)i?qi?p 

(C 8 ) 


+2cosv5pSinvJqSin(6»(p+,) - 9p)h^(_p_q)B^Bp^ . 


Appendix D: The in-in analysis for the full planar case 
In the following, we present the in-in formalism for the simplest case, which is the planar case. 


1. In-In integrals for the leading contributions in the full planar case 

In the following, by using the in-in formalism, we obtain the tensor-scalar-scalar bispectrum for each of the related 
diagrams, as given in the Fig. 


• SST: In this case, there is only one possibility, taking into account that we are using the first order in-in 
formalism here. Therefore, the TSS bispectrum would be. 


(^-(-kiCk2Ck3(??e)) = -24\/2B(^) sin^flfe, / ^ 


H 


'Vo 


/m(C(fc2, v)C{k 2 , ile)C{k 3 , r])C*{k 3 , 'qe)h+{ki,'q)h*_^_{ki,r]e)) + (^2 ^3^ 


/ V 2 RN, 
\ 


(^sin^9k, + P^^(ki)Pc(k 2 ) + (k 2 ^ ka)^ , (Dl) 


where 770 refers to an arbitrary time while rj^ denotes the end of inflation time. 

• SSV: Here we need to go to the second order in the in-in formalism. There are two different places for the 
interaction Hamiltonians to sit. They can either choose 771 or 772 . However, a detailed analysis shows that the 
final results would be the same for both of them so we can just calculate one of them and keep in mind that the 
final result needs an extra factor 2 , as has been written in the main text. 

/ jvf \ ^ /■''i / 

{h+kj^Ck^CksiVe)) = -96^21? I j sin^^fc^ / dru / dr]2i Im{h+{ki,r]i)h*^{ki,r]e)) 

\^Ve / J r)Q Jrjo \ 

Im(c{k2,m)C{k2,Ve)C{k3,7]2)C{k3,Ve)A'{ki,r]2)A *(^ 1 , 771 )^ -f {k2 ^ 

= - (3JV7 - 2JV.)sm“e.. ( n-(k2)P»Xki) 

+ (k2 ka)^ . (D2) 


• SVT: In this case, we need to consider the second order in-in formalism as well. Again there are two different 
locations for each of the Hamiltonian to sit. However, a detailed analysis shows that they will be the same. 
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Therefore in the following, we consider only one of them while keeping in mind that the hnal result would be 
twice, as is given in the main text 


Mf 


Hnl 


cm 


(/i+kiCk 2 Ck 3 (» 7 e)) = - 192 A/ 2 i? ( ) / dr]i j dr] 2 smek^[ smOk^^cosiOki - 0 k 2 )Im{C{k 2 ,m)C{k 2 ,Ve)) 

'Vo '^Vo 


Im(^C{h,r]2)Cik3,rie)h+iki,ri2)h*^{ki,r]e)A'{k2,r]2)A *ik2,r]i)'^ + (^2 fcs) 
( 2^/2R \ ^^^2 _ ( sin6»fci sin^fc^ cos (0^^ 


\ 


uz 

^3 


P^(k2)P?i+(ki) + (k2 <S4> ks)^ . 


(D3) 


• VVT: Here, one needs to expand the in-in formalism up to the third order in the perturbations. In this case, there 
are three different positions for each of the interaction Hamiltonians. However, as we have already mentioned 
above, these different locations lead to the same results, so we can only consider one of them here and remember 
that the final result needs an additional factor of three, as is given in the main text 

(^+kiCk2Ck3(??e)) = J J dri 2 J dr] 3 ^^^ ^sin0fc2 sin6»fe3 cos(0fcj - 

cos ( 6 »fci - 0 fc 3 )/m(C(fc 2 , Vi)Cik2,Ve))Im{C{k3,r]2)C{k3, r]e)) 
Irn[h+{ki,ri3)h\{ki,'qe)A'{k3,ri3)A *{k3,r]2)A'{k2,r]3)A *(A:2,??i)) + (^2 


- - Ne + sin sin 0^3 cos ( 6 »fc 3 - 0 ^ 3 ) cos {9ki - Oks) 

(^0Pc(k2)PH(ki) + (k2^k3)) . 


(D4) 


• SVV: Like in the previous case, one needs to consider the third order in-in formalism for this case. There are six 
different positions for the interaction Hamiltonians in this case. However, they would be the same and as in the 
above cases, it is enough to only consider one of them and keep in mind than the final result needs an additional 
factor 6 , as it is given in the main text 


(^+kiCk2Ck3(??e)) = - 384 V 2 P J dr]i J dr]2 J dr]3 ^ sin6*^3 sin6*^3 cos (0^^ - 6*^3) 

Im{h+{ki,r]i)h\{ki,T]e))lm{C{k2, m)Cik2,Ve)) 

Im(c{k3,T]3)C{k3,ile)A'{ki,r]3)A *{ki,'ni)A'{k2,'n3)A *(fc2,?72)) + (*2 k3)^ 

= - Sin Sin 0 fc 3 cos - 0 ^ 3 ) 

Pc(k 2 )P/^+(ki) + (k 2 4^k3)^ (D5) 


2. Proving the suppression of the terms without time derivative of the gauge field 

As it has been argued throughout the text, we do not consider terms which do not contain the time derivative of 
the gauge field. One intuitive reason might be due to the fact that magnetic field decays outside the horizon. So one 
can neglect them and only consider the electric held, as we did it for the whole of the analysis. However, in order 
to be more precise, it is worth to do a sample analysis for this case as well to have a better feeling about the order 
of magnitude of these terms. So in the following, we consider one example of the magnetic Held, for example let us 
consider the following interaction in the SVV sector 

L^agueticfield = V^ J d^P J d^qpqC{^p-q)ApAq (D6) 
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Using the in-in formalism, we have the following expression for the tensor-scalar-scalar bispectrum, 


(^+kiCk2Ck3(»?e)) = - (384v^i? 


/ Mp 


dm I dm j ^^3 ( — 1 ( fcifc2sin6ifej sinS 


\Hril 

/m(C(fc3, m)C{k3, Ve))Im{C{k 2 , m)C{k 2 , r?e)^'(fc 2 , m)^*{k 2 , m)) 
Im{h+{ki,ri3)hX{ki,'qe)A'(ki,m)A*{ki,m)) + (^2 k^) 

( 27V2r\ 


I lOOeH J 


sin0k^ sinOk^ (fci?7o) {k 2 'no) Pc0^2)Ph+0^i) + (k 2 <t4> ka) . 


(D7) 


Comparing Eq. (D7) with Eq. (D5), we see that the contribution from the magnetic held is indeed suppressed as 
compared to the electric held and it is safe to neglect it for the whole of our analysis. 


Appendix E: Detailed Analysis of the root-mean-square of the Quadrupole Asymmetry 

In this Appendix, we calculate the root-mean-square of the quadrupole anisotropy within our anisotropic setup. 
Starting from the dehnition of the root-mean-square of the Quadrupole asymmetry, 

= ( E iq2m|E (El) 

^ M^-2 ' 


we need to calculate the following quantity, 


I 2 

/l/f——9 \ / 


M=-2 

Now, in order to proceed, we can use the following formula for the GW 


/i*j(KL)/lah(KL)^ = P^(KL)M,,,b(KL) , 


(E3) 


where (Kl) is dehned in the following way, 

Mijabi^l,) = {PiaPjb + PibPja — PijPab) , (E4) 

in which Pij = (Sij — KpiKpj) is the projection tensor into the surface perpendicular to the GW wave number. 
Plugging back the above dehnition for the projection tensor inside Mijab{Kp), we would have. 


(^L)^ab(I^L)^ — Ph{,d^L){^iadjb dibdja dijS^^b ^ijKpb S^b^Li^Lj diaKpjKpb 

—SjbKpiKpa — dibKpjKpa — SjaKpiKpb + KpiKpjKpaKpb^ . (E5) 

We also recall the expression for the Aij(x,fi), 

Aij(x;, n) = — (n • x)Yn^) . (E6) 

_2 

Plugging the above expressions back into the Q , and using the following expression for the product of two YpM, 

YpM{n)Ylj^in')=(^^)pp{n-h'), (E7) 

M^-L ^ / 
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we would get, 

(f = {g^Nef Ph{KL) y d^n'P 2 (n • n')^l - (n • x)^ - (n'• x)^ + (n • x)(n'• x)(n • n') - 2 (Kl • x)^ 

+ (n' • x)(Kl • x)(Kl • fi') + (n • x)(Kl • x)(Kl • ri) - (n • x)(n' • x)(Kl • n)(KL • n') + (n' • x)^(Kl • x)^ 

+ (n • x)^(Kl • x)^ - (n • x)(n' • x)(n • n')(KL • x)^ + (Kl • x)'‘ - (n' • x)(Kl • x)^(Kl • n') 

-(n • x)(Kl • x)^(Kl • n) + (n • x)(n' • x)(Kl ■ x)^(Kl ■ n)(KL • n')^ • (E8) 

Now the remaining part is to calculate the above integrals. For this purpose, we use the following formulas, 


d^kPiik-h)Pi,{k-h') = 


) 6u'Pi{h ■ fi') 

2/ + 1/ ^ ^ 

1 + 1 


= (sTt)+( 2iTT) 


as well as Eq. (E7|, to rewrite Pi(n • ft') in terms of the products of YimS. 

We find the final result for the root-mean-square of the quadrupole anisotropy as. 


(E9) 

(ElO) 


q ' = {ST,){g,N,f Ph{KL)(l-{'Ki.-±f) 

/21iy/21l\/2 1 l\ 

l^OOOy q s J \ —t q' s' J 


(dOTT^) ^ ^ F*(x)n*,(x)F*,(KL)yiV(KL) 

= —1 t— — 2 


(Ell) 


As we see from Eq. (Ell), the final result is a function of Kl as well as the preferred direction. So we may also take 
the average of this function with respect to either Kl or the preferred direction, x. 
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